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Abstract 

We performed numerical simulations of 2D and 3D Edwards- Anderson spin glass 
models by using the recently developed multicanonical ensemble. Our ergodicity 
times increase with the lattice size approximately as V 3 . The energy, entropy and 
other physical quantities are easily calculable at all temperatures from a single sim- 
ulation. Their finite size scalings and the zero temperature limits are also explored. 

The low temperature phase of the spin glasses has distinct properties like broken 
ergodicity and the absence of self-averaging, which in turn make their numerical inves- 
tigation an extremely difficult task Q. For temperatures below a bifurcation point, the 
spin glass configuration space is supposed to split off into a number of valleys which are 
separated by high energy barriers. Due to the exponentially increasing relaxation times 
encountered in canonical simulations, tunneling between these many thermodynamics 
states becomes almost impossible. Recent simulations || of the 3D Edwards-Anderson 
model in a magnetic field seem to support the mean field picture @] rather than the 
alternative droplet model ||, but it can be argued that equilibrium at sufficiently low 
temperatures has not been reached ||. 

One of the simplest spin glass systems is described by the Edwards- Anderson Hamil- 
tonian 

H = - JijSiSj, (1) 

<ij> 

where the sum includes only the nearest neighbours and the exchange interactions = 
±1 between the spins s$ = ±1 are the quenched random variables. The constraint Jij 
is imposed for each realization. 

We present a new approach to spin glass simulations which reduces the exponential 
slowing to a power law and enables one to sample independent ground states in one 
simulation. This is achieved by exploiting the multicanonical ensemble 0. 
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The multicanonical ensemble can be defined by weight factors 



P M (E) = exp[-P(E)E + a(E)]. (2) 

a(E) and (3(E) are to be determined such that for the chosen energy range E m i n < E < 
-Emax the resulting multicanonical probability density is approximately flat: 

P mu (E) = c mu n(E)P M (E) « const (3) 

where n(E) is the spectral density In the present study we take E max = (/3(E) = 
for E > -E max ) and E min = E° the ground state energy of the considered spin glass 
realization. The purpose of the function a(E) is to give (3(E)~ l the interpretation of an 
effective temperature. This leads to the recursion relation 

a(E-A) = a(E) + [j3(E — 4) — /3(E)] E (4) 

where a(E ma , x ) = . The multicanonical function (3(E) is obtained via recursive multi- 
canonical calculations. One performs simulations with j3 n (E), n = 0, 1,2, which yield 
probability densities P n (E) with medians E^ edian . We start off with n = and /3°(E) = 0. 
The recursion from n to n + 1 reads 

/3 n (E) for E > E^ edian ; 



(3 n+l (E) = { 



(3 n (E) + 

0.25 x \n[P n (E + 4)/P n (E)\ (5) 
for E^ edian > E > E^ in ; 



(3 n+1 (E- iQ ) for E < E 



mm 



The recursion is stopped for m with E™^ 1 = E° being groundstate. 

By weighting with exp[— f3E + (3(E)E — a(E)\ canonical expectation values 0((3) = 
Z(/3)- 1 E E 0(E)n(E)exp(-(3E), where Z((3) = J2e n(E) exp(-@E) is the partition 
function, can be reconstructed for all j3 (the canonical inverse temperature). The normal- 
ization constant c mu in equation (3) follows from Z(0) = J2E n (E) = 2 N , where iV is the 
total number of spin variables. This gives the spectral density and allows to calculate the 
free energy as well as the entropy. 

First, we simulated the 2D Ising model by setting all the exchange interactions to 
Jij = +1 on a 50 x 50 lattice. Figure 1 stems from this simulation. Finite lattice energy 
density and specific heat results of Ferdinand and Fisher || were well reproduced. With 
Z(0) = 2 2500 used as input, we obtained 2.07 ±0.22 for the number of groundstates. After 
having this check, we considered the Edwards- Anderson spin glass. 

We performed multicanonical simulations of Edwards- Anderson spin glasses in 2D on 
L 2 lattices with L = 4,12,24,48 and in 3D on L 3 lattices with L = 4,6,8 and 12 0. 
We investigated as many as 32 realizations (having different configurations of quenched 
random variables Jy ) for smallest lattices and 4 realizations for the largest. Following 
an initial 2 x 10 5 sweeps, the thermal averages were evaluated over 4 x 10 5 to 8 x 10 6 
iterations depending on lattice size. 2 to 10 recursion steps of (5) were needed to reach 
to an approximately flat probability density (3). 

Our ergodicity time r£ is defined as the average number of sweeps needed to change 
the energy from £ max to E^ n and back. In Figure 2 the ergodicity time for 2D is plotted 



Figure 1: Ising model magnetic probability density. 



versus lattice size L on a log-log scale. A straight line fit gives r£ ~ £ 4 - 4 ( 3 ) sweeps, 
which in CPU time corresponds to a slowing down ~ V 3 ' 2( - 2 \ Here an exponential form 
yielded a completely unacceptable goodness of fit Q < 1CT 6 . 

The infinite volume groundstate energy and entropy are estimated from FSS fits of 
the form /j? = + c/V. For 2D these fits are depicted in Figure 3. Our groundstate 
energy density estimate e° = —1.394 ± 0.007 and groundstate entropy per spin estimate 
s° = 0.081 ±0.004 are all consistent with previous MC and the transfer matrix results 



TTfl. The reported groundstate entropy value translates into a large number of distinct 



groundstates, for instance, implies 1.1 x 10 81 groundstates on 48 2 lattice. 



Figure 2: Ergodicity times versus lattice size. 

For each realization we also evaluated spin glass susceptibility density Xq — (<? 2 ) /N 
and the Binder parameter B q = | 3 — (g 4 ) / (q 2 ) 2 . Our results were consistent with the 
previous evaluations [|TJ. 

For 4 2 lattice we calculated the exact results by enumerating all 2 16 configurations per 
realization and all quantities were found in agreement with the multicanonical evaluations. 

Next we compared our performance with multicanonical simulation to that of canoni- 
cal algorithm. For 12 2 and 24 2 lattices we performed an identical number of sweeps with 
multicanonical and canonical simulations for each realization. For all realizations the 
canonical tunneling rate is found smaller than the multicanonical one. When tunneling is 



Figure 3: FSS estimate of the infinite volume groundstate energy and entropy per spin. 



concerned, although it strongly depends on the realization, an estimated overall improve- 
ment factor 46 ± 22 for multicanonical simulations was observed for 12 2 lattices. For 24 2 
lattices we obtained zero tunneling events for most of the realizations in canonical simula- 
tions and about 10 - 20 tunnelings for the multicanonical. In view of the ergodicity times 
we achieved, this is no surprise. It clearly shows the superiority of the multicanonical 
approach. 

Dimension is a parameter in our codes. No difficulties were encountered in our in- 
vestigation of 3D spin glasses and results were encouraging ||. The groundstate energy 
density e° = —1.786 ± 0.003 and entropy per spin s° = 0.046 ± 0.002 are consistent with 
previous estimates ||12|| . For 3D we had a slowing down ~ ^/ 3 - 4 ( 2 ). 



For each realization we simulate two replica which differ only in the disordered starting 
configurations of spin variables Sj. Spin glass order parameter q can be defined as the 
overlap of the two replica || 

1 N 

m = mT,sIsi (5) 

To visualize the low temperature behaviour, we show in Figure 4 the spin glass order 
parameter distribution of one of the realizations on 8 3 lattice. 



Figure 4: Spin glass order parameter distribution on 8 3 lattice. 

Below the bifurcation temperature one clearly sees five configuration space valleys 
which are separated by high tunneling barriers. The multicanonical simulation overcomes 
these energy barriers by connecting back to the disordered high temperature states. The 
realization from which Figure 4 is depicted had nine tunnelings between the degenerate 



groundstates. For the Ising ferromagnet of Figure 1, the number of tunneling events was 
118. 

Our results make clear that the multicanonical spin glass simulations are very feasi- 
ble. Physical quantities can be evaluated at all temperatures from a single run. Tunneling 
through the high energy barriers and sampling of independent ground states become pos- 
sible, which certainly are relevant enrichments for the simulation of spin glasses. An 
investigation of the total weight and the overlap pattern of the spin glass order parameter 
might yield invaluable information about their groundstate structure [13]. 
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